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Abstract 

We analyze a model of interacting particles and strings described by a path integral with 
the Dirichlet boundary conditions. Such model is a natural framework to examine the pro- 
cesses involving the center-of-mass motion of string theory DO-branes: recoil, annihilation 
and pair production. We demonstrate that, within the proposed formalism, the exclusive 
annihilation/pair-production amplitudes admit a saddle point evaluation. Even though the 
saddle point equation cannot be solved analytically, it allows to extract valuable information 
on the coupling constant dependence of the amplitudes. In particular, DO-brane pair produc- 
tion turns out to be suppressed as exp[— 0(1/^)], much stronger than the naive expectation 
exp[— 0(1/ g st )}. All our derivations generalize rather immediately to the case of unstable 
DO-brane decay. In conclusion, we briefly comment on the possible implications our results 
may have for the conventional soliton-anti-soliton annihilation. 
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Introduction 



The physics of topological defects in relativistic field systems is a rather accomplished area of 
study PQ, yet fairly little remains known about how such localized objects are pair-produced 
in collisions of elementary particles, and how they annihilate. It is of very little use that 
the scattering of the fundamental fields quanta off a topological defect (a process related 
to annihilation by crossing symmetry) has been analyzed in quite some detail. The masses 
of the solitons are typically inversely proportional to the coupling constant, forcing the 
annihilation products into the high energy regions of the phase space. Transition between 
the scattering and annihilation kinematic regions would then require an analytic continuation 
non-perturbatively far in the phase space, and the perturbative expansion of the amplitude 
in the scattering region will be essentially unrelated to the values of the amplitude in the 
annihilation region. 

Even though there must exist a classical solution of the field equations describing the 
annihilation process, it does not appear to be accessible to any degree of analytic control. 
Indeed, there is apparently no straightforward way to prove that such a solution encapsulates 
even the most fundamental features to be expected from the annihilation process, e.g. the 
unitarity restrictions on the exclusive annihilation/pair-production amplitudes that we'll 
discuss below. 

One may expect that the situation becomes even more involved for the solitonic objects 
of string theory, the D-branes, since, in that case, we would have to deal with all the com- 
plexities of quantum gravity in addition to the problems that plague the treatment of the 
field-theoretical process. It has been known since the seminal work by Polchinski [2. that 
the solitons of string theory admit a particularly simple description in terms of the Dirichlet 
boundary conditions for open strings. Indeed, the analysis of elementary quanta scattering 
in the presence of D-branes is substantially simpler than the corresponding problem in the 
conventional relativistic field theory, as it merely amounts to computing the expectation 
values of certain operators in a free 2-dimensional field system. One may therefore hope 
that worldsheet techniques would provide a valuable tool for understanding the D-brane 
annihilation, a tool unavailable to the corresponding field theory considerations. 

Here, we attempt to construct a computational scheme for the DO-brane (D-particle) 
annihilation amplitude inspired by the simplicity of the CFT description of the extended 
D-branes. Since the center-of-mass motion of the D-particles is of a crucial importance for 
the processes we intend to consider, it is only natural that we introduce an explicit functional 
integration over the D-particle worldlines, to which the string worldsheets are attached via 
the Dirichlet boundary conditions. A similar method has been employed by Hirano and 
Kazama in their treatment of D-particle recoil [3]. It has been shown that, for low energy 
gravitational quanta, the scattering amplitudes obtained via this method can be matched 
with the standard no-recoil computation and respect the consistency requirements, such as 
BRST-invariance. Our goal in these notes is to address the question of how instructive this 
kind of approach can be for essentially high-energy processes, in particular, the D-particle 
annihilation. Similarly to Hirano and Kazama's paper, we restrict ourselves to bosonic D- 
particles, keeping in mind that the generalization to the potentially physically interesting 
case of superstring theory is likely to be conceptually straightforward, if only technically 
more challenging. 
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It is worth a notice that the attitude taken here is in many ways complementary to the 
attempts to address the same physical problem within the tachyon condensation paradigm 
of string field theory Indeed, the tachyon condensation approach typically deals with the 
initial state of coincident D-branes, whereas our aim is to analyze the genuine scattering 1 . 
On the other hand, should there be a need to generalize the approach advocated here to 
the case of higher dimensional D-branes, those will have to be wrapped around cycles in the 
compactification manifold, a requirement that does not appear within the string field theory 
considerations. 

Speaking of the possible phenomenological applications, one should certainly mention 
the kinetics of primordial topological defects. Even though the inflationary scenario suggests 
that the density of primordial topological defects is dramatically diluted during the period 
of exponential expansion [H], they could play an important role at the earlier stages of the 
life of the Universe 2 . The problem we're considering is also in a (somewhat indirect) relation 
to the microscopic black hole pair production, a subject of much theoretical controversy in 
the recent past 

Determining the dependence of the annihilation amplitude on the coupling constant is 
the central theme of our present investigation. It is a requirement imposed by unitarity 
[3 EHj that the exclusive annihilation amplitudes for the topological defects must be non- 
perturbatively suppressed, i.e. the amplitude should vanish identically when expanded in 
powers of the coupling constant 3 . The general intuition about how such suppression could 
arise dynamically appeals to the well-known property that the inverse mass of the topolog- 
ical defects is typically much smaller (by powers of the coupling constant) than their size. 
Therefore, the typical wavelength of the annihilation products will be much smaller than the 
size of the objects they are produced by, resulting in exponential damping. 

The problem with this kind of explanations is that the path integral includes arbitrar- 
ily singular configurations, and those could easily upset the exponential damping, which is 
supposed to underlie the non-perturbative suppression of the annihilation amplitude. (The 
precise meaning of this statement will become more apparent once we proceed with con- 
structing the actual formalism.) The saddle point estimate of the amplitude we obtain in 
the following sections can be seen as a proof that the measure of those singular configura- 
tions in the path integral is too small to alter the coupling constant dependence necessary 
to maintain unitarity. 

It is a fairly general principle that the processes involving topological defects upset the de- 
creasing significance of the final states with large multiplicity familiar from the perturbative 
Feynman calculus. Producing a large number of quanta may often turn out to be advanta- 
geous compared to the low multiplicity final states. A possibly more familiar manifestation 
of this kind of behavior is the emission of a large number of gluons in the instanton-induced 
processes of gauge theories fTT) . In our context, the bias for producing 0(1/ g st ) soft quanta 

*A scattering-like process in a tachyon condensation setting has been recently considered in 0. 

2 In the cosmological context, the medium in which the topological defects move may exert substantial 
influence upon the annihilation process 

3 Note that the total annihilation cross-section does not need to be non-perturbatively suppressed, and it 
is expected to be comparable to the square of the size of the topological defect for the case of gauge theory. 
This is one of the indications that the final state multiplicity distribution for topological defect annihilation 
is typically rather non-trivial. 
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is even more dramatic, as we'll be able to see after completing the saddle point evaluation 
of the amplitude. 

Quantization of D-particle worldlines 

When dealing with extended D-branes, one does not have to worry much about the center- 
of-mass motion: the infinite D-brane mass makes recoil impossible, at least when a localized 
object scatters off the extended D-brane. The situation is different, however, for compact- 
ified D-branes or genuine D-particles: for the low energy scattering experiments, the recoil 
is considerably suppressed by the mass of D-particle (inversely proportional to the coupling 
constant), but the center-of-mass motion becomes absolutely crucial in any essentially rela- 
tivistic process, such as annihilation. 

In the context of field theory, this problem has been addressed in ^2j. One introduces 
the translational (center-of-mass) modes for the soliton explicitly, and performs quantization 
with the center-of-mass coordinate treated as a canonical variable. It has been a subject of 
some controversy how the center-of-mass degrees of freedom should be introduced for D0- 
branes [T3J E] . Hirano and Kazama's recipe jSj constitutes a specific proposal to this end. 

To account for the motion of the D-particle as a whole, one introduces its coordinates 
explicitly and integrates over all the possible worldlines /^(i), with t being the proper time. 
The boundaries of the string worldsheet are restricted to the D-particle worldline, and the 
emission of closed strings is described by insertions of the closed string vertex operators in the 
interior of the worldsheet. Our main object of interest is the amplitude for two D-particles 
starting off at the positions x% and x% to annihilate into m closed strings carrying momenta 

G(xi,x 2 \ h, ■ ■ ■ , km) = £ / [Df] VtVX 6 (Jf M (0) - / M (t(0))) 

Vck J diff 

x exp [S D (f) - S st (X)} J] {g st V a (k a )} 

a=l 

where Sd is the action for the D-particle to be discussed below, S st is the standard (conformal 
gauge) Polyakov action 

the integration with respect to extends over all the inequivalent (unrelated by diffeo- 
morphisms) curves starting at x\ and ending at x 2 , the boundary of the worldsheet is 
parametrized by 9, and t(6) describes how this boundary is mapped onto the D-particle 
worldline. The sum is over all the topologies of the worldsheets (not necessarily connected, 
but without any disconnected vacuum parts) and x is the Euler number. Vck is the con- 
formal Killing volume (the negative regularized value of ^Hj should be used for the disk). 
The fully integrated form of the vertex operators is implied. We work in the Euclidean 
space-time, keeping in mind a subsequent analytic continuation to the Minkowski signature. 
The integration over moduli of the worldsheet is suppressed for the rest of this paper, as it 
does not affect the qualitative results. The annihilation amplitude can be deduced from ((H) 
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by means of the standard reduction formula: 
(ki, • ■ ■ , fc m |pi,p 2 ) 

= lim (pl + M 2 ) (pl + M 2 ) f dx 1 dx 2 e ipiXl e ip2X2 G{x 1 ,x 2 \k 1 ,---,k m ) (2) 
pI,p%-*-m 2 v ' v ' J 

where M is the D-particle mass. 

Several general questions have to be addressed regarding the expression (JTJ). Firstly, the 
issue with the Weyl invariance appears to be rather subtle. On the physical grounds, one 
would believe that making the D-branes fully dynamical reinforces the consistency of the 
amplitudes, much in the same way as respecting the supergravity equations of motion makes 
the non-linear a-models consistent. We shall examine below how this works explicitly to 
the lowest order of the recoil perturbation theory. The case of annihilation/pair-production 
is considerably less computationally straightforward, and assessing the issues of the Weyl 
invariance there is likely to require some new ideas. 

Of course, whether or not the integration over the D-particle worldlines reinforces the 
consistency of the string amplitudes depends crucially on the choice of the D-particle world- 
line action. It appears to be a fairly general principle JHI that the value of the effective 
action for a background which couples to strings is given by (minus) the sum of all con- 
nected vacuum string graphs evaluated in this background. Thus, very much in the spirit of 

s D [f}= £ { -^fjvtvx5(x,(e)-Utm)^vhs st (x)] (3) 

connected x 

Again, the negative regularized value of the conformal Killing volume should be used for the 
disk [15J. The exponentiation of the action in the path integral can be seen as a result of 
summing up the disconnected graphs containing vacuum parts 17J. It can be shown that, 
for nearly straight worldlines, the above action is reduced to the naive point-particle result 
M J dt. For curved worldlines, (J3J) would take into account the back-reaction from the space- 
time fields excited by the accelerating D-particle. Some properties of this action will become 
more apparent as we proceed with the computation of the amplitude. 

Given the close relation between the integrand in (0) and the worldline action (jSJ), it is 
convenient to rewrite (JIJ in the following form: 

G( Xl ,x 2 \ h, ■ ■ ■ , k m ) = / [VfjVtVX 5 (X,(6) - f^B))) 

all VcK J dlff (4) 

m 

xexp [-S st {X)\l[{g st V a {k a )} 

a=l 

Here, the sum extends over all the string diagrams including arbitrary disconnected vacuum 
parts. There is no explicit action for the D-particle worldline, but it arises after resumming 
the contributions from all the disconnected vacuum parts, which exponentiates to restore 
the So of 0. C is a combinatorial factor that can be deduced from (j2J). 

Using the transformation properties of the vertex operators under the target space trans- 
lations, it is easy to see that 



G(xi, x 2 \ h) = exp 



G 



X\ — X 2 X\ — x 2 
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The first term here merely provides for the momentum conservation 5-function in the Fourier 
transform, and (J2J) can be rewritten as 



lim ^ (pi + M 2 ) (pi + M 2 ) y dx exp 



;(pi -£2)2 



X X 



(5) 



We shall work with this representation in our subsequent calculation of the amplitude. 



The Gaussian integration 

By a direct inspection of it is easy to see that the integration over X is Gaussian and can 
be performed exactly. We will thus be able to recast the formalism into a (0+l)-dimensional 
form. The Gaussian integration we have to perform is closely related to the derivations in 
and can be implemented by applying the formula: 



j VX5{X{d) -£{0)) exp 



d 2 a 



Aiiot 



V1V1 + UX 



exp 



W / J(a)D(a,a')J(a')d 2 ad 2 a' -i J £(e)d n D(6,a')J(a')d9d 2 a' (6) 

' i{6)d n d n ,D{e,e')t{e')dede' 



Here, D is the Dirichlet Green function of the Laplace operator AD(er, a') = —5(cr — a'), 
and d n denotes the normal derivative evaluated at the boundary (which is parametrized by 
9). It is convenient to consider 



2' 2 



J 



E 

all 



^£ / [Vf] VtVX8(X,(9) - /„(«(*))) 
Vck J diff 



x exp 



-S st (X)+i J d 2 aJ^ 



(7) 



instead of (jlj). Indeed, differentiating with respect to the source J and setting it to ki5{cr — 
<7j) allows us to reproduce the amplitude for an arbitrary vertex operator insertion. Perform- 
ing the integration in (J7J) by means of (jUJ) yields: 



G 



J 



E 

all 



C (9st) x 



V, 



CK 



exp -ixa J J»{a)D{a, a')J^a r )d 2 ad 2 a' 



I [vf] m m exp H / f^KWnDie^j^ded 



V 



(8) 



exp 



Ana' 



7 / r(t(e))d n d nl D(e,e')f ll (t(o , ))dOd6' 



It should be noted that D(a, a'), x, C and V depend on the topology of the diagram corre- 
sponding to each particular term in the sum. For diagrams with disconnected parts, D(a, a') 
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is block diagonal in the sense that it vanishes whenever the two arguments belong to different 
disconnected components. 

The path integral in (JHJ) may appear rather cumbersome, as one of the functions to be 
integrated over appears in the argument of the other one. Nevertheless, the integration over 
can be performed exactly by means of a technique very similar to the treatment of the 
free point-like particle in [T8] . 

We first rewrite the measure on reparametrization equivalence classes of f^{t) as 



Pf\m = Vf& f -1 =Vf Vzexp 



z(f 2 - l)dt 



where 5[f 2 — 1] is a product of 5-functions at every point (reinforcing t to be the proper time), 
and, for each t, the integration over z(t) is along a contour going from c — ioo to c + ioo 
in the complex z plane, with c being an arbitrary (positive) constant. This contour can 
of course be deformed, an opportunity implicit in our subsequent application of the saddle 
point method. 

If we now introduce 

Af(t, t') = J d9d9'd n d n ,D{e, 6') 5 (t- t{6)) 5 (f - t{6')) , d(t, a) = J d6d n D(6, a)5(t- t{6)) 
the /-integration in (jHJ) can be recast into a manifestly Gaussian form: 



2' 2 



j) = ^^M! exp - na > J r(a)D(a,a')J^a')d 2 ad 2 a' 



all VcK 
Vt Vz Vf exp 



- / z{f - l)dt 



cxp 



-i I f^(t)d(t,a)J f ,(a)dtd 2 a' 



exp 



4na' 



- / nt)Af(t,t')f,(t')dtdt' 



It is convenient to pass on to the integration over Vf by means of the relations 



Vf = dTVf 5 I J fdt + x» 

vO 



X 



F(t) = - + J Fdt 





We should also keep in mind that 

T 



J dW(t, t') = J dt d(t, a) = -1 
o o 



as 



J d6 d n d n 'D(6, 9') =0 J d6d n D(6,a f ) = -1 
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If we also perform the Fourier transform of (0), we arrive at the following representation 



G(pi,p 2 \J) = J dxexp 



G 



x x 
2'~2 



J 



all 



exp 



W / J»(a)D(a, a')J^a')d 2 ad 2 a' 



J dTVtVz el zdt J Vf exp i J f^t) (rf - j d(i,a)J»(a)did 2 a 



dt 



exp 



- F(t)B(t,t')f,(t')dtdt> 



where we've introduced 



t v 



B(t, t') = z(t) 5(t - t') + — j dt J dt'Af(t, i') 





At this point, the Gaussian integration becomes completely straightforward and yields 



G( PuP2 \J) = Yl 



all 



V c 



exp 



K 



-ixa' J J»{a)D(a, a')J^a')d 2 ad 2 a' 



dTVt(9)Vz(t) det[i3]- 1/2 exp 



z dt 



(9) 



exp 



t V 

\ I ^ ~ I c/ (^ CT ) JAt ( CT )^ 2cr ) i3_1 (^ t/ )(piA t - / d(i\a')J^a')dt'd 2 a')dtdt' 



From this representation, it is apparent that the endpoints of the z integration contour can 
be moved towards — oo. The contour itself cannot shrink to — oo, however, on account of the 
singularities of (detB)" 1 ^ 2 . Should there be a discontinuity in t(6), these singularities move 
towards — oo allowing the contour to be deformed arbitrarily far to the left in the complex z 
plane. Then the integrand will vanish due to the last factor in the second line of Q. This is 
as it should be, since discontinuous worldsheets do not give any contribution to the original 
path integral (HJ). 

The fact that we have been able to perform the integration over the D-particle worldlines 
exactly is rather remarkable, since, for any given worldline, the action Sd featured in (fT|) 
cannot be computed. This action however is itself constructed from string amplitudes and 
can be completely eliminated from the path integral at the cost of including the disconnected 
string worldsheets, as it has been done in (0J. This peculiar simplification plays an impor- 
tant role in making possible the saddle point evaluation of the annihilation/pair-production 
amplitudes. 



The saddle point method 

The presence of the p\ and J in the exponential of the last line of © alludes to the relevance 
of the saddle point techniques, as those quantities become non-perturbatively large in the 
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annihilation kinematic region. This can be made more apparent by introducing pi — Pi/M 
and J = J /M (M being the D-particle mass) and rewriting Q as 



g(pu P 2\j)=y: 



all 



C (9st) x 

V C K 



exp 



cxp 



M 2 



na'M 2 J J»(a)D(a, a')J^a')d 2 ad 2 a' 
dTVt(6)Vz(t) det[£]- 1/2 exp 
Pi- f rf(t»J rM (ff)^V)i3- 1 (t,t , )(pi / , - / d(t' \a')J^a')dt'd 2 a')dtdt' 



z dt 
t 



o 



The latter representation is still somewhat inconvenient, since the position of the saddle 
point depends on the saddle point parameter M. This dependence is very simple, however, 
and can be completely eliminated by reseating 4 z{t) — > Mz(t/4:7iMa'), t{6) — ► 4nMa't(9), 
T — > 4irMa'T, which brings the above integral to the form 



St 



G(pi,p 2 \J) = J2 7} ' ex P 

VCK 



all 



-na'M 2 J J"(a)D(a, a')J^o')d 2 od 2 o' 



dTVt(d)Vz(t) det[£]- 1/2 exp 



exp 



4WM 2 / : dt 
t' 



(10) 



vra'M 2 J (pi- J d- J»did 2 a)A- l (t, t') (p lfl - j d- J^dt' d 2 a')dtdt' 
o o 



A(t, t') = z(t) 8(t -t')+ / dt / dt' Af(t, t') 



where 



In view of the subsequent application of the saddle point method, it is necessary to specify 
how the integrand is analytically continued to the complex t(8). It may seem worrisome 
that the definitions of Af(t, t') and d(t, a) involve ^-functions, but these 5-functions appear 
in integral convolution, and the result can be made manifestly analytic by expressing them 
through their Fourier series on the interval [0,T]. Such a representation defines Af(t,t') and 
d(t, a) as periodic with respect to shifting t or t' by multiples of T. It is worth a notice 
that analogous analytic continuation is not possible at the level of (jHJ) due to the lack of 
smoothness 5 in / M (t). In this respect, the summation over all the worldlines performed in 
the previous section is a crucial prerequisite for the saddle point techniques to be applicable. 

With these specifications, the saddle point equations can be readily derived by requiring 
the sum of the two exponents in (|10|) to be stationary with respect to variations of t(6) and 

4 It should be kept in mind that, from now on, the geometrical parameters of the worldsheet in physical 
space are related to those represented by t(9) through a factor of 0(l/g s t)- 

5 It is a familiar fact (related to the central limit theorem) that the measure on random paths is dominated 
by curves of fractal dimension 2: see, for example, |19j . 
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z(t). In this manner, one obtains 



t(6>) 



v»{t{e)) Ude'd n d n ,D(e,e') J v tl (t')dt' + J d n D(e,a)j^(a)d 2 a \ =o (11) 



where we've introduced 



V^Vy, = -4 



/' / t' 



(12) 



V»(t) = J dt'A-\t,t') \Pi- J d(t \a')J»(cj')dt'd 2 G 
\ 

Multiplying (fTT|) by 5(t — t{9)) and integrating with respect to 9, we get 



zV 2 + -z 







which, coupled with ()12|). implies z = const. The relevant value of the constant can be deter- 
mined from the following argument. Because of the subsequent application of the reduction 
formula, the value of the amplitude is determined by the divergent piece of G(pi,p 2 \J) as the 
momenta go on-shell. This divergence can only come from large values of T in the integral 
(which, of course, corresponds to the long distance propagation of the on-shell particles). 
For large values of T, there will always be parts of the worldline arbitrarily remote from the 
location of the worldsheet. There, (fT2|) will fix z = ±i\fp I /2 (with <p 2 = p\ = pi). Since the 
saddle point value of z does not depend on t, we conclude that, for the purposes of evaluating 
the amplitude, z can be set to ±iy/p^/2 everywhere. 

Thus, we have two complex conjugate saddle points for z. After substituting the cor- 
responding values to (fTTj) . we should obtain two complex conjugate solutions for t{9). The 
equation (fTTj) is clearly beyond the reach of analytic methods, but the saddle point configu- 
ration of the disconnected vacuum components of the worldsheet can be easily identified as 
t(9) = const. Indeed, for 9 belonging to a disconnected vacuum component, the second term 
in the brackets of (fTTj) vanishes. Enforcing t(9) = const makes the first term in the brackets 
vanish as well, the equation being thereby satisfied. 

Moreover, the values of the second functional derivative of the saddle point functional 
with respect to t(9) can be calculated for the parts of the worldsheet which do not emit any 
final state particles. Indeed, if 9 and 9' both belong to a disconnected vacuum component 
with t{9) = t , 



5 2 



5t(9)5t(9>) 



1 1 
J (pi+ J d-J^dtd^A- 1 ^^')^^ J d-J^dt'dV 



-2V 2 {t )d n d n ,D{9,9') = 8d n d n ,D(9,9') 



Furthermore, if 9 belongs to a disconnected vacuum component and 9' belongs to any other 
component of the worldsheet, the corresponding second derivative vanishes. 
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The above specifications in fact allow to perform a complete resummation (to the leading 
order of the saddle point approximation) of the contributions to (|1U|) coming from the discon- 
nected vacuum components. Indeed, for each disconnected component in a given diagram, 
we'll obtain a factor 



The value of this integral is a constant 6 (depending on the topology of each particular 
disconnected vacuum component) times a factor of T, which comes from the integration 
over the constant mode of t(9). If we now recall that the combinatorial coefficients C in (jTUj) 
originated from expanding the exponential in (0), it is easy to see that, to the leading order 
of the saddle point approximation, the effect of all the disconnected vacuum components 
amounts to a factor of exp[/iT], where \i is a constant, which can be evaluated by computing 
the Gaussian integral (fT3j) . In fact, it is easier to notice that G(pi, pi\J = 0) is merely the 
D-particle propagator, and, identifying its pole with the D-particle mass M, conclude that 
H = -A-koi'M 2 . 

For the parts of the worldsheet which do emit final state strings, it is not possible to 
solve the equation (jllj) explicitly. However, valuable information can be extracted from the 
mere assumption that a saddle point exists (this can be checked in the low energy scattering 
case). It is an important circumstance that, for the worldsheets located not too close to 
the endpoints of the worldline, in the limit T — > oo, the value of the saddle point function 
does not change under the shifts of t{6) by a constant (we assume p\ = pi = p 2 , as it is 
on-shell). It is not hard to take into account this quasi- zero mode though, as the subsequent 
application of the reduction formula discards everything but the term growing most rapidly 
as T goes to oo. Since extending the interval to which t and t' belong to [—00, 00] does not 
introduce any singularities to (JTTJ), we should expect that, as T goes to 00, the solutions to 
(fTTj) located not too close to the endpoints of the worldline approach a fixed shape. The 
saddle point evaluation of the Dt and Dz integrals in (fTUj) can be therefore written (up to 
the pre-exponential factors) as 



where T does not depend on T and can be determined from the solution to the equation (fTTj) 
with t and t' running from —00 to 00. The above expression is written for the saddle point 
z = i\/p*~/2, the contribution from z = —iy/p I /2 is the complex conjugate of this. We shall 
omit the o(T) term in the following, as it vanishes upon the application of the reduction 
formula. 

Assembling everything together, up to the pre-exponential factors, the saddle point esti- 

6 One may argue that, for higher genus worldsheets, the modular integrations (that we do not write 
explicitly) will exhibit a tachyonic divergence. The same problem arises if one tries to compute the higher 
genus corrections to the D-particle mass. Both problems will disappear in the superstring case, which we see 
as a conceptual rationale behind the formal identification of the value of the integral with the D-particle 
mass, which we're about to make. 




(13) 




(14) 
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mate of (fTUjl is 

G( Pl ,p 2 \J) 



(gstf 



saddle Vck 



exp 



-na'M 2 J J"{a)D{a, a')J^a')d 2 ad 2 a' 



x 



4toV 1 + ii/p 2 T 



<iT + c.c. 



where the summation is once again performed only over those worldsheets which do not 
contain any disconnected vacuum parts. Performing the T integration and applying the 
reduction formula, we finally obtain 



x (exp {-a'M^fJ/M]} + c.c. 



(15) 



To compute the actual amplitude, one has to combine this expression and its functional 
derivatives with respect to J so as to imitate the vertex operator insertions, as well as to 
integrate over the worldsheet moduli and the positions of the vertex operators (these will be 
saddle point integrations themselves, in analogy to [201 )■ All these tasks are likely to require 
numerical computations, as is the evaluation of the ^"-functional. 

An essential piece of information contained in the expression (|15|) is, however, that it 
identifies the dependence of the annihilation amplitude on the coupling constant. We can 
see that, for a final state with the number of quanta fixed as g s t — > and momenta growing 
proportionally to the mass of the D-particle, the amplitude diminishes as exp(— 0{\/g 2 st )). 
Indeed, the anticipated non-perturbative suppression has become manifest 7 . It is worth 
noting that the character of this coupling dependence is not at all obvious at the intermediate 
stages of our computation. Only upon the application of the saddle point method can we 
see that the qualitative arguments outlined in the introduction do, in fact, result in a non- 
perturbative suppression of the amplitude. 

It may appear rather surprising that the estimate of the amplitude obtained above is 
exactly the same as if we had used the simplistic free particle action M J dt in place of the 
expression © which takes into account the emission and re-absorption of virtual string states 
by the accelerating D-particles. The qualitative explanation here is that, being objects of 
effective size v^o 7 , the D-particles are not likely to emit virtual string states of energy much 
higher than l/y/a'. These states will produce worldline curvatures of order g st and will not 
therefore affect the leading order of the saddle point approximation. One may have suspected 
that the worldline curvature could become large in the region where the final state strings 
are emitted. Yet, the saddle point equation (jllj) implies that the amplitude is dominated by 
worldsheets of physical size ~ l/g s t (which corresponds to t(9) being of order 1). Therefore, 
the emission occurs from a segment of the worldline whose size is of order l/g s t, and the 
relevant curvatures are of order g st again. 

7 There doesn't appear to be any straightforward way to prove that T is greater than zero in the entire 
annihilation kinematic region (after the analytic continuation to Minkowski signature). Nevertheless, (|15|) 
suggests that the amplitude is either non-perturbatively suppressed or non-perturbatively enhanced, and the 
inherent absurdity of the latter option is quite apparent. 
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Recoil perturbation theory and the Fischler-Susskind 
mechanism 



Let us now take a step back and examine how the more conventional coupling expansion 
of the recoil perturbation theory arises in our set-up. Namely, let's set aside the (neces- 
sarily relativistic) annihilation/pair-production processes and consider closed string states 
scattering off a D-particle. Since the mass of the D-particle diverges in the limit g s t — > 0, it 
can be treated as static to the lowest order in g st (if we keep the momenta of the incident 
closed strings fixed as g st — > 0), and the corrections due to the motion of the D-particle's 
center-of-mass (i.e. recoil) will appear as a perturbative expansion in powers of g st . This is 
the familiar recoil perturbation theory. 

The formalism of the previous section is, of course, perfectly applicable in this case. Since 
the momenta of the closed string states (and hence J) are kept fixed as g st — > (and hence 
M — > oo), J = J/M can be treated as perturbation. The saddle point equation (|TT|) implies 
then that the saddle point configuration t(6) = const + 0(1/M), i.e. the path integral (fTU|) 
is dominated by nearly constant t{9) (the fluctuations of t(9) are of order 1/M generically, 
and so is the saddle point value, as we've just remarked). One could look for solutions 
of the saddle point equation as an expansion in J . However, it is much more efficient 
to generate the perturbation theory directly from the path integral ([TUjl by constructing a 
suitable expansion of the "effective action" functional of ()10|) 



t f 
+WM 2 J (fi - J d- J^didtyA'^t,?)^ - J d- J^di' d 2 a')dtdt' 
o o 

There is one difficulty one encounters in implementing such a program. It is a most 
straightforward approach to try to construct a Taylor-like expansion of S e ff in powers of 
t{6) around t(9) = const: 



S eff [z(t),t(e)] = -Ana'M 2 / zdt 



S e ff (const + t(9)) 

= S e// (const) + Jddf^ t{e) + l _ J d0d0 > * S «f 



5t(6)5t(6') 



t{9)t{6') + 



This strategy comes to mind in immediate relation to the computational techniques most 
commonly used in cx-models, and it has been employed in [3 for the purposes we're presently 
pursuing here. Unfortunately, such an expansion does not exist. In £3 . various ^-function 
prescriptions have been devised (in the lowest order) to deal with the infinities arising when 
one tries to brute-force the Taylor-like expansion, but the situation quickly becomes hopeless, 
if one tries to envisage the general structure of the recoil perturbation theory in such a 
framework. 

The origin of the above complication can be traced back to the non-analytic properties 
of the worldlines in the path integral (jHJ. Indeed, for non- analytic f^it) (most worldlines 
are fractal and therefore non-differentiable jTHj) the "effective action" in (JSJ) does not admit 
a Taylor-like expansion in t(8) around any configuration of t(9). The integration over f^{t) 
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improves the situation considerably: the resulting effective action can be expanded around 
any t{9) 7^ const, but the non-analyticity still survives for the worldsheets whose boundary 
shrinks to a single point. (One may become worried about whether such non-analyticity could 
undermine the application of the saddle point method from the previous section. Whereas 
any discussions of explicit contour deformation in multi-dimensional case are necessarily 
rather subtle, there are still good chances that the deformation implicit in the saddle point 
evaluation of the integral can be attained, since S e ff is analytic for any t(9) 7^ const.) 

Luckily, the Taylor-like expansion is not the only way to generate a sensible perturbation 
theory. Appearing as insertions in a Gaussian path integral, the exponentials of t(9) are just 
as tractable as powers of t(6). We shall therefore resort to a combination of a Taylor-like 
and a Fourier-like expansion. We first introduce 

A(t, t') = i^-bit - t') B{t, f) = 5z(t) 5{t -t') + J dtj dt'Af(t, ?) 







such that A(t,t') = A(t,t') + B(t,t'), and expand formally 

(we work with one of the two saddle points in z(t), the other one will give a complex conjugate 
contribution, as in the previous section). Upon inserting this expression in the S e ff we can 
isolate the following terms that need to be retained in the exponent (everything else is small 
and may be treated perturbatively, as it will become apparent after we exhibit the powers 
of the expansion parameter 1/M): 

S^ f = Ama'M 2 ^T 
t t' 

Sf) f = -Ana'M 2 J dtdt' J dt J dt'Af(t,t') = -Ana'M 2 J t{0)d n d n ,D(6,9')t{9')d9d6' 







gfj , brio* [ 
VP 



t 



S±% = 47U ^fV [dt [d- J,dtd 2 a = ^^-'A I t(8)dnD(0, a)J{a),ddd 2 a 
v P J J V P J 

(17) 

The first term here merely provides for the correct pole structure. It will disappear after 
integration over T and application of the reduction formula. The remaining three terms 
define a Gaussian integral with respect to z(t) and t(9). 

A general term in the expansion (fTBj) can be constructed as follows. We should take a 
certain number of factors 

U U+i 

J dt J dt'N{t,t') (18) 


and, for each i add a factor 

(zVP72)-"* +1 
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where rii are non-negative integers. Then, for the first and last t% we either add a factor of 
pi or a factor of 

J d{t,a)Ji"{a)dtd 2 a (19) 
o 

We then integrate over all ij's and multiply the result by ira'M 2 . 

We now proceed with the Gaussian integral over 8z(t). As it can be seen from the expres- 
sion for Sg A, each factor of 5z (t) will produce a factor of 1/M in the result of the Gaussian 
integration. Where several 5z's occur at the same ti, divergencies from the singularity of 
the 5z propagator (proportional to the (^-function) will be present. Just like in the case of 
free particle considered in [IB] , these divergencies will merely renormalize the saddle point 
(expectation) value of z(t). 

We are left with a Gaussian integral over t(6). As we've remarked before, the structures 
of the type ([TB]) and ([T9*]l cannot be expanded in powers of t{6) around t{6) = const. Instead, 
we'll Fourier-trasform in t,'s: 



o o 

t 



dt f dt'Af{t,t') -> / e^^dnd^Die^'y^+^dede' 

U 

f d{t,a)J»{o)did 2 (j -f - / e iK > m d n D(9,a)J»(a)d9d 2 a 
o 

At this point, the integration over t(8) can be explicitly performed. The substitution t{6) — > 
t(6)/M, Ki — > Mki reveals the additional powers of 1/M that each term will receive upon 
the integration over t(9). 

Having sketched the general structure of the recoil perturbation theory, we shall now 
turn to the important issue of how the consistency of the string amplitudes is respected 
to the lowest non-trivial order. As we've remarked early on in the course of our present 
investigation, the inclusion of curved D-particle worldlines into the path integral (0) will 
generally induce a Weyl anomaly thereby threatening the decoupling of the negative norm 
states, which is a crucial requirement for the consistency of the string theory S-matrix. The 
educated hope is that (a suitable modification of) the familiar Fischler-Susskind mechanism 
|21[ will come to rescue the formalism that has been constructed so far. Indeed, it is a 
well known fact that there are singularities in the modular integration for the higher genus 
worldsheets coming from the corners of the moduli space. Should we have chosen the D- 
particle worldline action correctly, the Weyl anomalies induced by those singularities will 
precisely cancel the Weyl anomalies coming from the coupling to an accelerating D-particle. 
Let us see how this works out to the first non-trivial order of the recoil perturbation theory. 

The Fischler-Susskind mechanism in our set-up implies a cancellation between an infrared 
divergence on a higher genus worldsheet and an ultraviolet divergence on a lower genus 
worldsheet (in the sense of cancellation of the corresponding Weyl anomalies). To the lowest 
order, this means that we should examine the ultraviolet divergencies on a disk to the order 
1/M and compare them with the modular integration divergencies from an annulus coupled 
to a straight D-particle worldline (the curved worldlines will only contribute to the higher 
orders in 1/M). 
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When we expand S e f / (for a disk) according to ()16|) and treat all the terms except for the 
ones indicated in (|17|) as a perturbation, we can identify the following insertions that can in 
principle contribute to the order 1/M in the resulting recoil perturbation theory Gaussian 
integral: 

2ma'M 2 



8WM 2 



Vp 2 



dt / d(t,a)J^(a)dtd 2 a 



t t' 



\ dtdt' / dt / dt'Af{t,t') / dt"d 2 ad{t",a)J^{a) 



o o 



1 (8na'M 2 



^ \ dtSz(t) / d(t,a)J tl (a)dtd 2 a 



P 

After we integrate over 5z(t) and Fourier-transform, these three terms become: 



, I dn 
Ft ~ a I — 



e lKt(e) d n D(9,a)J^a)d9d 2 a 



a'Mpf 



dn 



d9d9't(9)d n d' n D(9,9')e iKt ^ / e~ iKt( - e "^d n D(9" ,a)J^a)d9"d 2 a 



dn 



Pi 



F3 ~ a 

We now have to compute the path integral 

Vt(9) (F 1 + F 2 + F 3 ) exp 



$ f e iKm d n D(9,a)J IM (a)d9d 2 a 



o(2) r.( 4 ) ' 

The relevant structure that is featured in all the three terms is 
dn 



iK{t{e)-t{e')) 



x exp 



Vt{9) e 



-Ava'M 2 / t(9)d n d n ,D{9,9')t{9')d9d9' + 



47ua'M 



t(9)d n D(9,a)J^(a)d9d 2 a 



a'M 



exp 



WpiPi / J tl {a)d n D{9, a)[d n d n/ D]- 1 {9, 9')d n >D{9', o')J v {p')d9d9' d 2 od 2 o 

f dn 
x / -^-exp 



-2k 2 loge + k/ii(^, 9') + k 2 /i 2 (#, 0'; 

Where loge = d n d n /D(9,9) is the regularized value of the singularity of the boundary-to- 
boundary propagator (e being the worldsheet cut-off), and hi t 2(9,9') are certain (cut-off 
independent) functions. If we now substitute k = n/^loge, we observe that the functions 
hi t 2 do not contribute to the UV divergent piece of the path integral 8 , which turns out to be 
(9, #')-independent and proportional to 

\/log£ 



a'M 



iViPi / J M (ff)a n D(e,ff)[a n a n ,i)]-He,e')o n /-D(e', ( T')J l ,((T')<iede' ( i 2 CT dV 



(20) 



8 The IR divergence at k = is not what interests us here, and it is in fact unphysical. It must disappear 
in a more accurate treatment taking into account that t only varies within the finite range [0, T]. 
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Of course, the k integral above is divergent for small finite values of e, and the expression 
we've given must be understood in terms of analytic continuation in e. As we'll see below, 
an analogous analytic continuation is needed for the annulus modular integration. This is 
very similar to what has been described in |14j . 

Now, in the path integral of F 2 , the expression (|2()jl will appear in convolution with 
d n d n 'D(9,9'), and, since the latter is orthogonal to constant modes, the path integral of F 2 
will not contain any UV-divergent piece. Both path integrals of F\ and F 3 will simplify due 
to the relation 

/ d n D(6, a)J^(a)d6d 2 a = - J J»d 2 a = Vi ~ V2 

and the UV divergent part of the path integral of F3 will turn out to be of order 1/M 2 , not 
1/M due to the relation 2pi(pi — p 2 ) = (pi — p 2 ) 2 ■ If we also notice that the exponential 
in (120)1 is merely the value of the Gaussian integral 



Vt(6) exp 



Q(2) q(i) 

'^eff <=>eff 



(without any insertions), we conclude that the UV divergent part of the disk amplitude to 
the order 1/M is proportional to 

(Pi -P2) 2 



where A Q D is the zeroth order value of the disk amplitude, i.e. the amplitude computed in 
the background of a straight D-particle worldline, in neglect of recoil. 

We now have to compare this result with the divergence in the modular integration of the 
annulus amplitude to the lowest order in 1/M, i.e. the annulus amplitude in the background 
of a straight D-particle worldline. The computation of such an annulus amplitude is fairly 
standard. The fact that the disk amplitude A° D makes an appearance in the expression 
(|2*T|) should be taken optimistically as far as the prospective cancellation of divergencies is 
concerned, since the familiar plumbing fixture construction of j22] relates the divergencies 
in the modular integration to amplitudes on worldsheets of lower genus. Let us see in some 
more detail how this actually works out. 

Along the lines of [22| 5 the annulus amplitude with an insertion of the operators V^ l \ • ■ ■ , V^ n ' 
(in the interior) can be expressed through the disk amplitudes with additional operator in- 
sertions at the boundary as follows: 



(yW ■ ■ ■ V^) annulug = £ / dqq h "- 2 J dOdO' (V a (6)V a (e')V^ . . . 



where the summation extends over a compete set of local operators V a {6) with conformal 
weights h a , and q is the gluing parameter that can be related to the annular modulus. The 
divergence in the integral over q coming from the region q « (corresponding to an annulus 
degenerating into a ring) will be dominated by the terms with the smallest possible h a . 
Neglecting the tachyon divergence, which is a pathology peculiar to the case of the bosonic 
string, we identify (in close relation to the investigations of JH]) the following set of relevant 
operators (parametrized by an integer i 6 [1, 26] and a real to): 



V\e,u) =:d n X\6)e 



iujx (e) 
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with conformal weights h = 1 + a'uo 2 /A. For small values of q (which is the region we're 
interested in) only small values of to will contribute into the integral. We therefore transform 
the annular divergence as follows: 



where we've taken into account that the operator / d n X(9)d6 merely shifts the (straight) 
D-particle worldline and inserting it into any amplitude amounts to multiplication by the 
total transferred momentum p\ — p 2 . If we now cut off the dq cicu-integral at q — e, the 
regularized value is proportional to ^/\oge. Since q is introduced as a distance along the 
boundary of the worldsheet, the geometrical meaning of e here is precisely the same as in 
the formula (J21|) . Recalling that there is an additional power of g st in the annulus amplitude 
as compared to the disk amplitude with the same insertions, we conclude that the modular 
integration divergence is proportional to 



To the order 1/M (i.e. to the order g st ), we recognize precisely the same structure, as 
in (|21jl. We shall not work out here the actual coefficients that depend on a number of 
conventions, but it is clear that the non-trivial dependencies on the momenta and the cut-off 
do in fact match, as it is necessary for the successful implementation of the Fischler-Susskind 
mechanism. 

One cannot help but noticing that the version of the Fischler-Susskind mechanism that 
we've just described bears a strong resemblance to the investigations of [Hj. Indeed, limiting 
themselves to the lowest order of the string coupling expansion, the authors of ^1] have shown 
that, to that particular order, the consistency of the string S- matrix should be restored if the 
background of a straight D-particle worldline is augmented by an inclusion of the operator 



(where Q(X°) is a step function). Heuristically, such an operator corresponds to the D- 
particle abruptly starting to move with the appropriate recoil velocity at the moment X° = 0. 
(Once again, the operators / d n X l (9)d6 with i being one of the Dirichlet directions shift the 
entire trajectory of the D-particle.) The UV divergence induced by this operator is precisely 
what we've found in ()21j1 . Moreover, the divergent integral 



makes an actual appearance in the derivations of [T^j. These parallels should not be too 
surprising since, for example, the structure 




9 st (Pi ~ Pi) 



2 (y« . . . vM) = g st (p a - p 2 f A^fi^e 






o 
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(entering the expression for the saddle point value of the D-particle trajectory develops 
a discontinuity in its first derivative as the boundary of the worldsheet shrinks to a single 
point. 

The authors of [Hj have remarked that the abruptness in the change of the D-particle's 
velocity inherent to their scenario should probably be alleviated at the higher orders of the 
recoil perturbation theory. In fact our present formalism bypasses this problem from the very 
beginning. The summation over all possible D-particle trajectories automatically gives rise 
to the UV divergence required by the Fischler-Susskind mechanism, making it completely 
unnecessary to introduce the abrupt recoil by hand. 

Is it possible to reach a deeper understanding of how the Fischler-Susskind mechanism 
works in our present formalism beyond the lowest order of the recoil perturbation theory? 
Extending our above analysis to the higher orders in g st is likely to require a more system- 
atic picture of both the general structure of the recoil perturbation theory and the general 
structure of the divergencies in the modular integration. Yet, conceptually, there doesn't 
seem to be anything puzzling about how the cancellations could actually occur. 

Of course, the situation is much more intricate for the annihilation/pair-production case. 
In that regime, the first non-trivial UV divergence will arise if we try to compute the pre- 
exponential factor in the saddle point estimate (|15j). There doesn't appear to exist any 
meaningful expansion of this pre-exponential factor in any small parameter. It is merely a 
determinant of some integral operator, whose divergent part needs to be matched against the 
modular divergencies from the higher genus worldsheets. Any cancellations in this setting 
will necessarily be rather subtle. In particular, they may require certain relations between the 
values of the saddle point exponentials, and therefore certain relations between the solutions 
to the saddle point equation (fTTj) for worldsheets of different genera. We shall not pursue 
this line of thought any further in this present investigation. Even though the considerations 
related to the Fischler-Susskind mechanism are of a crucial importance for establishing the 
consistency of the worldline formalism, they do not affect the estimates of the saddle point 
exponential and, in particular, its coupling constant dependence, which is the main result of 
this paper. 

Emission of a large number of strings 

One may be somewhat worried by the fact that the annihilation amplitudes computed via the 
saddle point method are suppressed much stronger in the limit g st — > than exp(— 0(1/ g st )), 
which is generally expected to be the order of the non-perturbative corrections in string 
theory (221- The resolution of this apparent paradox is that the actual annihilation is most 
likely to proceed into a large (~ l/fl'st) number of quanta, whereas the number of strings 
in the final state has been kept fixed as g s t — > in our preceding considerations. Indeed, 
as we could see in the above, the annihilation process was dominated by worldsheets of size 
~ l/g s t, and it would hardly be natural for such extended configurations to decay into a 
small number of strings, a circumstance responsible for the anomalously strong suppression 
of the amplitude (|I5J). 

We shall now obtain some estimates of the amplitude for the D-particles to annihilate 
into a large number of quanta, and show that those contributions are in fact likely to saturate 
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the bound of [2o]. The analysis of the previous sections still applies, but the position of the 
saddle point (and hence the functional J 7 ) now depends on g st , and therefore extra care needs 
to be taken in making statements about the coupling constant dependence of the amplitude. 

It can be seen by inspecting the expression (JTKJ) that, unless all but an infinitesimal 
fraction of the final state strings are emitted from disconnected worldsheet components 
of disk topology with exactly one vertex operator on each of them, the g* t factor itself 
will provide suppression of order exp(0(ln(g st )/ g st )). The extent to which the dynamics 
(accounted for by the saddle point estimates) enhances this suppression depends on the 
topology of each particular diagram. When most of the final state strings are emitted from 
disks with a single vertex operator, the non-perturbative suppression can only arise from the 
value of the saddle point exponential. We shall now estimate the contribution from those 
diagrams. 

When the worldsheet consists of many disconnected components each of which carries 
away only a small fraction of the momentum, the independent displacements of those dis- 
connected components cause only a small change of the saddle point function. We therefore 
encounter quasi-zero modes very similar to the ones described above (jl4j) . One way to ac- 
count for these quasi-zero modes is to introduce the integration over them explicitly into 
the path integral by inserting S(J i t(9)d6 — U ^d6) for each of the disconnected components 
of the worldsheet (the 9 integral runs over the ith component of the worldsheet, U thereby 
being the average position of its boundary). The variational problem of (|!Hll2jl has then to 
be solved subject to the constraints specified by the 5-functions. (This would introduce the 
corresponding Lagrange multipliers into the right-hand side of ([lip, but they are inessential 
to our subsequent considerations.) 

In order to obtain constraints on the value of the amplitude, let us expand the V(t) in 
(jlljl in Taylor series around fj (for each i). It is easy to see that unless more than 0(l/glf £ ) 
(with e arbitrarily small and positive) of the t^s come into an interval of size 0(g st ), the 
solution to (fTTj) for 9 belonging to the ith disconnected component is given by 9 

t(9) =U + \J [dnd n >D{9, 9')}- 1 d n >D(9', a')V^t i )J lM {a')d 2 a'd9' + 0(g l + £ ) 

(The second term here is closely related to the saddle point which dominates the recoil 
perturbation theory in the lowest order of g st .) For the value of the saddle point function 
(the exponential in the last line of (|10p). we obtain 

/ T t(6) t{8') 



-na'M 2 



J dtz{t)V 2 {t) + J d9d9'd n d n ,D{9,9') J V{t)dt J V{t')dt' 

\0 

2nia'M 2 J^T 



ira' 



Y, I d9d9 , d 2 ad 2 a'V^i l )J fl (a)d n D(9,a) [d n d n 'D(9, 9')]" 1 d n /D(9', a')V v (Ji)J v {a') 



+0(l/glr" 



9 We are essentially demanding the variation of V to be small over the extent of a single disconnected 
component of the worldsheet, and then use the expression for the saddle point value of t{6) in the background 
of a straight worldline. 
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The second line here merely provides for the correct pole structure (and disappears after the 
application of the reduction formula). Each term in the sum from the third line is of order 
0{9st) an d negative (V is imaginary to the lowest order in g st ). However, since there are 
0(1/ g s t) terms, we should expect a non-perturbative suppression of order exp(— 0(l/g st )). 

We still have to deal with the excluded region of integration over the tj's, namely, the 
configurations with more than 0(1/ g\t £ ) of the Vs in an interval of size 0(g st ). The volume 
of this region is proportional to exp(0(hn(g st ) / gli £ )) , however, and this is the amount of non- 
perturbative suppression provided regardlessly of the value of the saddle point exponential. 
Barring any unfathomable cancellations, the sum of the contributions from the two regions 
would be of order exp(— 0(l/g s t)), in accordance with the general estimates of the non- 
perturbative effects in string theory. 

Bosonic D-particle decay 

It is a notable feature of the worldline formalism that most of the results obtained in the 
previous sections are immediately generalized to the case of bosonic D-particle decay that 
has received a large amount of attention within the tachyon condensation considerations. 
Indeed, the natural proposal for the D-particle decay amplitude is 

(fa, ■ ■ • , k m \p) — lim / dxdx' e ipx G(x,x'\ fa, ■ ■ ■ , k m ) (22) 

p 2 ^— M 2 J 

where G(x,x') is as it has been defined in (0). The meaning of this formula is that the 
endpoint of the D-particle worldline is unconstrained in the path integral, and the reduction 
formula is applied only to the starting point. The path integral thus describes a "disappear- 
ing" D-particle. 

Given the formal similarities between the expressions (j2J and (J22J), it is not surprising 
that the derivations presented above for the annihilation/pair-production case carry over to 
the decay process with a formal substitution p 1 = p, p 2 = 0. One important difference is 
that, for p\ ^ p\ the quasi-zero mode described in the paragraph above (fTlj) will be absent. 
The saddle point configuration of the worldsheet will be localized near the endpoint of the 
D-particle worldline, as one would expect on general grounds. For that reason, the factor 
of T featured in (jl4j) will not make an appearance for the case of D-particle decay. This 
is reassuring, since it precisely replaces the double pole necessary for the application of the 
reduction formula (J2J) by a single pole necessary for the application of the reduction formula 
(J22|) . The conclusions regarding the coupling constant dependence of the amplitude to the 
leading order of the saddle point approximation will remain intact, i.e. the exclusive decay 
amplitudes will be non-perturbatively suppressed as exp[— 0(l/g 2 st \. 

It would be interesting to see to which extent the results obtained here are compatible 
with the investigations of the closed string emission by a decaying DO-brane within the 
tachyon condensation approach 4J. Unfortunately, the (non-self-consistent) neglect of the 
closed string back reaction that underlies the derivations of |3] makes it hard to judge which 
features of the emission spectrum described there should be trusted. 
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Conclusions and speculations 



As the calculations of the previous sections reveal, the problem of DO-brane annihilation 
appears in many ways much more tractable in our present context than the field-theoretical 
soliton-anti-soliton annihilation or, more generally, any other topological defect annihilation 
process that the author is aware of. Indeed, the saddle point method provides rigorous 
results regarding the coupling constant dependence for the annihilation with a fixed number 
of final state quanta, and certain qualitative estimates can be devised for final states whose 
multiplicity increases as g s t goes to zero. For a fixed number of final state closed strings (and 
for DO-brane pair production), the amplitudes turn out to be suppressed as exp[— 0(1/ g^ t )], 
much stronger than the naive expectation exp[—0(l/g st )}. 

For the case of closed string scattering off a D-particle, the worldline techniques developed 
here offer a rather elegant construction of the expansion in string coupling constant that 
properly takes into account the D-particle recoil. In particular, we resolve the paradox 
related to the abruptness in the change of the D-particle velocity that made an unwelcome 
appearance in the previous attempts to address the same problem |T3j . 

It is a somewhat dissatisfying feature of our present formalism that only exclusive am- 
plitudes appear to be computationally accessible. Should one have had actual D-particles 
at one's disposal, exploring the refined properties of their annihilation (such as final state 
multiplicity distribution) would certainly provide a valuable insight into their nature. Exclu- 
sive annihilation amplitudes would be the relevant theoretical prediction to guide this sort of 
experimentation. However, in most phenomenological applications (such as cosmology), the 
most important quantity is the simplest characteristic of the annihilation process, namely, 
the total cross-section. 

It is the non-trivial nature of the final state multiplicity distribution, i.e. the dominant 
role of the final states with a huge number of quanta, that makes the reconstruction of the 
total cross-section from the exclusive annihilation amplitudes so hard. If one could overcome 
this difficulty, it would be possible, for example, to examine in our setting the stretched D- 
string interconnection, a process that has so far only admitted a semi-quantitative treatment 
|24j . The generalization of our present formalism to the case of D-strings is not likely to 
pose much difficulty. The problem is that, in the cosmological context, one is interested in 
the total interconnection probability. 

For the case of D-particle pair production, we have established a very strong non- 
perturbative suppression (of order exp(— l/g^/)). It would be interesting to contemplate 
whether it may have anything to do with the exponential suppression of the microscopic 
black hole pair production amplitudes argued by one of the sides in the dispute jS]. 

Another amusing feature is that there does not appear to be a unitarity relation that 
constrains the absolute normalization of the D-particle annihilation amplitude. It is hard to 
judge at this point whether or not this represents an actual non-perturbative ambiguity in 
string theory. 

Since many qualitative features of the topological defect annihilation, such as the non- 
perturbative suppression of the exclusive amplitudes and the high mean multiplicity of the 
final state, appear to be rather general, it is natural to ask whether our computation can 
shed any light on the field-theoretical soliton annihilation, a process that has thus far defied 
an analytic treatment. The Christ and Lee formalism of ^2] provides a rather conspicuous 
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link to our approach, as it introduces explicitly the center-of-mass position of the soliton, a 
worldline to be integrated over. It is certainly true that performing the path integral over the 
fluctuations of the fields around the given trajectory of the soliton is a much more difficult 
task than the worldsheet integrations pertinent to the D-particle case. This circumstance is 
also likely to undermine the direct applicability of the saddle point method, since it seems to 
hinge upon the summation over the worldlines. Yet, one could try to develop some intuition 
about the kinds of worldlines that contribute most, and establish some bounds, say, on the 
coupling constant dependence. If such a program succeeds, it would effectively provide a 
path-integral derivation of certain non-trivial properties of the classical annihilation solution. 
Mathematically and physically, such a derivation would present a fair amount of elegance. 
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